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Abstract. The systematic investigation of ground-state shape evolution from 
7— unstable 0(6) to spherical U(5) for even— even 112 - 134 Te has been presented by 
using the quadrupole moment constrained Hartree-Fock-Bogoliubov (HFB) method 
with the Skyrme force SLy4. 124 Te has been pointed out as to be the possible critical- 
point nucleus with E(5) symmetry. 
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1. Introduction 

Spherical vibrator, rotational ellipsoid, and deformed shapes of nuclei are related with 
modes of collective motion and geometric shapes of nuclei [H [21 El S] - When the num- 
ber of protons and neutrons is changed, nuclei can show changes of their energy levels 
and electromagnetic transition rates among collective modes of nuclei. Transitions from 
one kind of collective behavior to another are named shape phase transition. These 
transitions are quantum phase transitions (QPTs) [5]. They are different from thermal 
phase transitions which occur as a function of temperature. This implies changes in 
the shape of the nucleus. The control parameter is the number of nucleons. In the last 
decade, many researchers have given insights into the evolution of structure of nuclei 
(in particular transitional regions of rapid change) They have used the concepts 
of QPTs, phase coexistence and the critical-point symmetries (CPS) [TJ IE] proposed by 
Iachello, as well as a raft of geometrical models [91 [101 HH [121 (131 EH]- 

Theoretically, QPTs have been mostly studied in the Interacting Boson Model 
(IBM). It holds the U(5), SU(3) and 0(6) symmetries within the simplest U(6) symme- 
try. The U(5), SU(3) and 0(6) dynamical symmetries correspond to the shape phase of a 
spheroid, axially prolate rotor and 7— soft rotor, respectively [15]. Using the model, the 
authors pointed out that second order shape phase transition occurs between U(5) and 
0(6) while first order shape phase transition occurs between 0(6) and SU(3) [HI ITT] . 
More recently, Iachello has pointed out that the properties of nuclei lying at the critical- 
point of a shape phase transition can be described by special solutions of the Bohr 
Hamiltonian, called as critical-point symmetries. The E(5) critical-point symmetry [TJ 
corresponds to the second order critical-point between U(5) and 0(6), while the X(5) 
critical-point symmetry [8 J corresponds to the first order transition between U(5) and 
SU(3). Experimentally, The E(5) and X(5) symmetry has been realized in the spectrum 
of 134 Ba [18] and 152 Sm [19], respectively. The introduction of critical-point symmetries 
E(5) and X(5), and their experimental realizations have triggered many works on quan- 
tum phase transitions [20J. 

In spite of the fact that the IBM and the solutions of Bohr Hamiltonian are use- 
ful for determining shape phase transitions in nuclei, mean field formalisms (e.g., HFB 
method [2j EEl 22] and relativistic mean field (RMF) model [231 EH ESI ES]) which pro- 
vide a correct prediction of many nuclear phenomena have been successfully used to 
study the shape phase transition in nuclei. The RMF theory has been used to inves- 
tigate the critical-point nuclei in even— even Sm [2T] and Ce [28] isotopes. In these 
studies, 148 . 150 . 152 Sm and 128 . 13 °. 132 > 134 Ce have been suggested as an example of the pos- 
sible critical-point nuclei with X(5) symmetry. Beside, Ti isotopes have been examined 
in the HFB method [29] and RMF model [30] to investigate the critical-point nuclei. 
In these studies, 48 > 52 > 60 Ti and 46 > 52 > 60 Ti have been found as to be the possible candi- 
date critical-point nuclei with E(5) symmetry in the RMF model and HFB method, 
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respectively. A series of isotopes in rare-earth region have been found as to be pos- 
sible critical-point nuclei [3TJ [32]. Also, Mo isotopes were investigated by using the 
RMF theory [331 El] and 94 Mo has been suggested as to be 7— unstable nucleus in 
Ref . [33] • In these studies, potential energy curves (PECs) come by quadrupole moment 
constrained calculations have been examined in order to identifmg of the critical-point 
nuclei. Relatively flat PECs correspond to the critical-point nuclei with E(5) symmetry, 
while PECs with a bump correspond to the X(5) symmetry (The relation between shape 
phase transition and PECs can be found in [HI [35]). However, for a quantitative analy- 
sis of QPTs in nuclei electromagnetic transition rates and ratios of excitation energies 
should be calculated [36]. For this reason, the generator coordinate method (GCM) has 
been used to perform configuration mixing of angular-momentum and particle-number 
projected relativistic wave functions restricted to axial symmetry in [37J. In recent 
years, the GCM has been extended on triaxial states [381 ESI HQ]. Nevertheles, the ap- 
plication of these methods in a systematic study of shape transition is at present still 
very time-consuming because of its triaxiality. It should be noted, however, that the 
potential energy curves obtained from quadrupole moment constrained calculations are 
important, and can provide a qualitative understanding of the shape phase transition. 
Particularly, the evolution of the PECs along the isotopic or isotonic chains can be use- 
ful for investigation of shape phase transitions in nuclei. 

Rapid shape changes in nuclei have been known for about half a century [6]. Clas- 
sic shape transition regions take place in the light Si-Mg region [H], near A = 100 
(Z ~ 40) [12], light rare-earth region (A ~ 150) and actinides. Beside, the 7— unstable 
character of nuclei in the mass region A = 120 — 130 was pointed out in [43] many 
years ago. And, 124 Te in this region was experimentally investigated and suggested to 
be a possible 7— soft nucleus [44J. In this study, constrained HFB method has been 
used to obtain the ground-state properties of even— even 112 ~ 134 Te isotopes such as total 
binding energy and quadrupole deformation. The shape evolution of Te isotopes has 
been analized by examining their PECs. 



2. The theoretical framework and numerical details 

In the HFB formalism, a two-body Hamiltonian of a system of fermions by means of a 
set of annihilation and creation operators (c, c^) is given by 

n\U2 n 1 n2nyrn 

where v nin , 2n . ini = (nin 2 \ V \ n^n^ — n 4 n 3 ) are anti-symmetrized matrix elements of 
the two-body N-N interaction. The ground-state wave function |$) is described as the 
quasi-particle vacuum a>k\&) = and the linear Bogoliubov transformation: 

a k = Y,{U*kCn + V* k cl), al = Y,(VnkC n + U nk cl) (2) 
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provides connection between the quasiparticle operators (a, cv) and the original particle 
operators. The basic building blocks of the theory are the density matrix and the pairing 
tensor. In terms of the normal p and pairing k one-body density matrices: 

/W = ($k,c n |$> = (V*V T ) nn ,, n nn , = ($\c nlCn \$) = (V*U T ) nn , , (3) 

the expectation value of the Hamiltonian ([ID could be expressed in terms of an energy 
functional: 

($\H\$) 1 1 

E[p, k] = = Tr[(e + -I>] - — Tt[Ak*] (4) 

where r niri3 S«2^4 ^nxn^n^n^prun^ and Z\ nin2 ^^nzn^^nxn^n^n^n^n^- Modern energy 
functional (T4]) includes terms that cannot be simply related to some prescribed effective 
interaction |45j . In terms of Skyrme forces, the HFB energy (Tjf]) has the form of local 
energy density functional: 

E[p,p] = jd 3 H(r), (5) 

where H(r) = H(r) + H(r) is the sum of the mean field and pairing energy densities. 
The variation of the energy ([S]) according to the particle local density p and pairing 
local density p results in Skyrme HFB equations: 



E 



/ h(r, a, a') h(r,a,a') W U(E, ra') 
v h(r,a,a') -h(r,a,a') ) [ V(E,ra') 

E + X W U{E,ra) \ 
E-X) J { V(E,ra) J ' 



(6) 



where A is chemical potential. Local fields h(r, a, a') and h(r, a, a') can be calculated in 
the coordinate space (2j S6] . 

The HFB equations flH]) have been solved by expanding quasi-particle wave functions 
that conserve axial symmetry and parity on a harmonic oscillator basis expressed in 
coordinate space prescribed by Stoitsov et al. jl6] in the present study. For pairing, 
the zero-range pairing interaction is taken into account and Lipkin-Nogami method 
is employed (Further details are given in Ref. [46J). The oscillator parameter 6 is 



chosen as to be bo = y2(/?, 2 /2m)(49.2y4 _1 / 3 ). In order to obtain the potential energy 
curves (PECs) in the present study, the standard quadratic form of the quadrupole 
constraint [211 HS] has been performed. The standard quadratic form can be interpreted 
by the formula E Q = Cq((Q) — Q) 2 , where Cq is the stiffness constant, (Q) is the 
average value of the mass quadrupole moment operator (Q = 2z 2 — r 2 ) and Q is the 
constraint value of the quadrupole moment. For describing the deformation of nuclei, the 
quadrupole deformation parameter is commonly used rather than quadrupole moments. 
The relation between the constraint quadrupole moment Q and constraint quadrupole 



deformation parameter f3 2 is given by the formula f3 2 = \Jtt/5(Q)/ (r 2 ) jl6]. 16 oscillator 
shells have been taken into account in the present calculations. There can be found 
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a number of effective Skyrme forces in literature for correct prediction of the nuclear 
ground-state properties of nuclei [47J HH H9] . In this work, the Skyrme force SLy4 [49] 
has been used to calculate ground-state properties of even— even Te isotopes with 60 
< iV < 82. 

3. Results and discussions 

The calculated total binding energies for even— even 112 - 134 Te isotopes obtained from 
the constrained HFB method with the SLy4 parameter set are tabulated in Table [TJ 
Also, the predictions of RMF model [50] and the experimental data [51] are listed for 
comparison. The total binding energies of all isotopes are reproduced well by the SLy4 
Skyrme force. The deviations are at most 0.3%. Also, as can be seen in Table [TJ the 
predictions of RMF model with NL3 interaction are in a good agreement with the ex- 
perimental data. The mean differences between experimental data and the predictions 
of the HFB method and RMF model with NL3 interaction are 2.532 and 1.688 MeV, 
respectively. For this reason, it can be pointed out that the predictions of both HFB 
method and RMF model are good in describing the ground-state binding energies of 

112-134 rp e 

The mean field formalism based on the Hartree-Fock approximation with phe- 
nomenological effective interactionsis is important in the microscopic description of 
nuclei [521 [53]. It allows a unified description of the ground-state properties for nu- 
clei throughout the nucleidic chart. One of the great achievement of the theory is that 
not only it can reproduces binding energies and densities, but it also provides a good 
description of the size of the ground-state deformations in nuclei [54]. In this study, 
the quadrupole deformation parameters 02 for 112_134 Te have been obtained from con- 
strained HFB method with SLy4 Skyrme force. They are shown in Fig. [TJ Also, the 
predictions of RMF model with NL3 interaction [50] and the experimental data [55] are 
shown for comparison. As can be seen in the Fig. [TJ the calculated f3 2 values obtained 
from HFB method with SLy4 parameters are in good agreement with the experimental 
data. Only amplitude of quadrupole deformation parameter 2 values obtained from 
both of the HFB method and RMF model are given in the Fig. [TJ and the exact values 
of 02 are listed in Table [2] It should be noted, however, that 2 cannot be observed 
directly from an experiment. To obtain the 02 value from an experiment, a conventional 
way is that the electric quadrupole transition rate from ground-state + to the 2 + state 
B(E2) t can be used [55J. The correlation between B(E2) f and 2 is given by the 
formula 2 = (4tt/3ZR%)[B(E2) f /e 2 ] 1/2 where R = 1.2A 1/3 . The formula based on 
rigid rotor cannot always represent a parameter of deformation. The extracting of 2 is 
questionable in the case of spherical nuclei, because B(E2) f connects vibrational states 
in the spherical nuclei. In particular, the radius Rq is so small for light nuclei. This 
elicites a very large 02 deformation with the formula. However, in medium-mass and 
heavy region usage of the formula is siutable [30J. 
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In this study, the constrained HFB method with Skyrme force SLy4 is employed 
for the investigation of shape evolution of even— even 112 - 134 Te because of its success 
in describing the binding energies and quadrupole deformation parameters for Te iso- 
topes. In Fig. El the potential energy curves (PECs) for 112_134 Te. In the figure, the 
total binding energy of Te isotopes for the ground-state is taken as to be the reference. 
In the Fig. [21 starting from 112 Te to 120 Te, the nuclei have oblate shape. In the PECs 
of the 122 ~ 126 Te, their barriers against deformation are weak which means that these 
nuclei may be in a transitional region. In particular, the PEC of 124 Te in the Fig. [2] 
seems flat from f$i = —0.2 to 02 — 0.25. Through these 02 ranges, the variation of 
the binding energies in the PEC of 124 Te are less than 0.4 MeV. This implies that the 
barriers against deformation are very weak, and 12 Te maybe a possible example of the 
critical-point nuclei with E(5) symmetry. With increasing of neutron numbers start- 
ing from the 128 Te to 132 Te, shape of Te isotopes have become prolate and finally 134 Te 
which has shell closure with magic neutron numbers N = 82 are found to be as spherical. 

In Table El the differences of the binding energy between the spherical-state and 
the ground-state of even-even 112 ~ 134 Te isotopes are shown to understanding of how the 
shape of the Te isotopes changes with the neutron number as an additional evidence to 
the results of the PECs. They can show how soft the nucleus is against deformation. The 
calculated binding energy differences between the spherical-state and the ground-state 
of 112_134 Te isotopes changes from to 2.572 MeV. Drastic changes are clearly visible 
around 120_124 Te in the binding energy differences. There is a clear jump appearing 
at 124 Te which implies that 124 Te can be a possible candidate critical-point nuclei with 
E(5) symmetry. 

As an additional evidence for confirmation of the result of this study, the 
ratios of experimental excitation energies of 124 Te nucleus |44j are given in Table HJ 
Also, the U(5), X(5), SU(3), E(5) and 0(6) symmetry predictions are listed for 
comparison [71 El Hi]- The characteristic ratio -R 4 / 2 = E(4f)/E(2f) and the ratio of the 
energies of the first two excited + states R0/2 = E(02)/E(2f) are tabulated. As can 
be seen in Table IH the E(5) symmetry values obtained from solution of Bohr-Mottelson 
differential equations for R A / 2 and i? /2 is 2.20 and 3.03, respectively. They are closer 
to the observed ratios i2 4 / 2 = 2.07 and R / 2 = 2.75 which means that 124 Te may hold 
the E(5) symmetry. 



4. Summary 

The total binding energies and quadrupole deformation parameters for even— even 
ii2-i34rfe isotopes have been calculated in the constrained HFB method with Skyrme 
SLy4 force as in a good agreement with experimental data. The ground-state shape 
evolution of these nuclei are investigated by using the potential energy curves. 124 Te 
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has been found as to be an example for possible critical-point nucleus, which marks 
the phase transition between spherical U(5) and 7— unstable shapes 0(6). It should be 
noted, however, that this work investigates only the PECs of Te isotopes with respect 
to the P degree of freedom. For a quantitative identifying of the E(5) symmetry in Te 
isotopes, one should go betond mean field. 
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Table 1. The total binding energy for the ground-state of 112 134 Te in units of MeV. 



This work RMF [5U] Exp [5T 



112 Te 


937.821 


938.880 


940.610 


ii4 Te 


958.058 


959.060 


961.337 


ii6 Te 


977.630 


978.530 


980.860 


118 Te 


996.680 


997.600 


999.454 


120 Te 


1014.148 


1015.640 


1017.281 


122 Te 


1030.893 


1032.530 


1034.333 


124 Te 


1047.472 


1049.160 


1050.686 


126 Te 


1063.296 


1066.980 


1066.368 


128 Te 


1078.858 


1080.750 


1081.439 


130 Te 


1094.204 


1096.430 


1095.941 


132 Te 


1108.850 


1112.220 


1109.914 


134 Te 


1123.508 


1126.430 


1123.435 
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Table 2. The ground-state quadrupole deformation parameter /3 2 for Te isotopes. 



This work RMF [50] Exp 



112 Te 


0.187 


0.164 




ii4 Te 


0.304 


0.232 




iie Te 


-0.171 


0.257 




118 Te 


-0.173 


0.175 




120 Te 


-0.169 


0.179 


0.201 


122 Te 


-0.135 


0.161 


0.185 


124 Te 


-0.096 


0.138 


0.170 


126 Te 


-0.093 


-0.003 


0.153 


128 Te 


0.076 


-0.002 


0.136 


130 Te 


0.062 


0.032 


0.118 


132 Te 


0.028 


0.000 




134 Te 


-0.005 


0.000 
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Table 3. The difference of the total binding energy (in units of MeV) between the 
spherical state and the ground-state of 112 - 134 Te obtained by the constrained HFB 
method with SLy4 Skyrmc force. 



lNFiiplpi 

1 i I. IA. ±v~± 


HFB-SLv4 


112 Tp 

-Lt: 


9 ^9Q 


114 T 




116 T 


9 ^1Q 


118 Te 


2.572 


120 Te 


1.900 


122 Te 


1.259 


124 Te 


0.354 


126 Te 


0.410 


128 Te 


0.227 


130 Te 


0.160 


132 Te 


0.053 


134 Te 


0.000 
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Table 4. The ratios of available experimental excitation energies for 124 Te isotopes 
with some theoretical predictions for comparison [6—8, 44]. 





R4/2 


Ro/2 


U(5) 


2.00 


2.00 


X(5) 


2.91 


5.67 


SU(3) 


3.33 


>2 


E(5) 


2.20 


3.03 


0(6) 


2.50 


4.50 


Exp 


2.07 


2.75 
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Figure 1. The ground-state quadrupole deformation parameter for Te isotopes. The 
predictions of the HFB method with SLy4 Skyrme force are compared with the those 
of the RMF model with NL3 interaction and experimental results. 




Figure 2. The potential energy curves for 112 134 Te obtained from the constrained 
HFB method with SLy4 Skyrme force. 



